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Abstract 

The general and explicit relation between phase times and dwell times for quantum tunneling 
or scattering is investigated. Considering two identical propagating wave packets symmetrically 
impinging a one-dimensional barrier, here we demonstrate that these two distinct transit time 
definitions give connected results where, for such a colliding configuration, the phase time (group 
delay) accurately describes the exact position of the scattered particles. The analytical difficulties 
that arise when the stationary phase method is employed for obtaining the phase (traversal) times 
are all overcome since the multiple wave packet decomposition allows us to recover the exact 
position of the reflected and transmitted waves. In addition to the exact relation between the 
phase time and the dwell time, leads to right interpretation for both of them. 
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Finding a definitive answer for the time spent by particle to penetrate a classically for- 
bidden region delimited by a potential barrier has occupied the physicists for decades [1-12] . 
In particular, people have tried to introduce quantities that have dimension of time and can 
somehow be associated with the passage of the particle through the barrier or, strictly speak- 
ing, with the definition of the tunneling time. These proposals have led to the introduction 
of several transit time definitions that can be summarized by three groups. (1) The first 
group comprises a time-dependent description in terms of wave packets where some features 
of an incident packet and the comparable features of the transmitted packet are utilized to 
describe a quantifiable delay as a tunneling time [13]. (2) In the second group the tunneling 
times are computed based on averages over a set of kinematical paths, whose distribution 
is supposed to describe the particle motion inside a barrier. In this case, Fcynman paths 
are used like real paths to calculate an average tunneling time with the weighting function 
exp [i S x{t)/h], where S is the action associated with the path x{t) (where x{t) represents 
the Feynman paths initiated from a point on the left of the barrier and ending at another 
point on the right of it [14]). The Wigner distribution paths [6], and the Bohm approach 
[15, 16] are included in this group. (3) In the third group we notice the introduction of 
a new degree of freedom, constituting a physical clock for the measurements of tunneling 
times. This group comprises the methods with a Larmor clock [3] or an oscillating barrier 
[17]. Separately, standing on itself is the dwell time defined by the interval during which 
the incident fiux has to exist and act, to provide the expected accumulated particle storage, 
inside the barrier [18]. In spite of no general agreement [12, 19] among the above definitions, 
the so called phase time [20] (group delay) and the dwell time have an apparently well estab- 
lished relation between them [13, 21]. However, these time definitions remain controversial 
since in the opaque barrier limit they predict effective tunneling velocities that exceed the 
vacuum speed of light and may even become unlimited (Hartman effect) [22]. For instance, 
some of the barrier traversal time definitions lead, under tunneling conditions, to very short 
times, which can even become negative. It can precipitately induces an interpretation of 
violation of simple concepts of causality. Otherwise, negative speeds do not seem to create 
problems with causality, since they were predicted both within special relativity and within 
quantum mechanics [8]. A possible explanation for such time advancements can come, in 
any case, from consideration of the very rapid spreading of the initial and transmitted wave 
packets for large momentum distribution widths. Due to the similarities between tunneling 
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(quantum) packets and evanescent (classical) waves, exactly the same phenomena are to be 
expected in the case of classical barriers [36]. The existence of such negative times is pre- 
dicted by relativity itself based on its ordinary postulates [12], and they appear to have been 
experimentally detected in many works [23, 24]. In some recent analysis, the proportional- 
ity between the phase time and the time averaged stored energy was used to explain these 
unusual effects [11, 21]. In particular, from the time-independent Schroedinger equation, a 
relation between the group delay and the dwell time was derived for quantum tunneling . 

In this manuscript we are concentrated on overcoming the difficulties that appear when 
the stationary phase method (SPM) is utilized for deriving tunneling times in order to 
accurately investigate the relation between the phase time delay and the dwell time for 
quantum tunneling or scattering. Taking into account the restrictive conditions for the use 
of the method, we report about a theoretical configuration involving a symmetrical collision 
between two identical wave packets and a one-dimensional rectangular potential barrier 
[25]. Using the procedure we call multiple peak decomposition [26], we demonstrate that 
summing the amplitudes of the reflected and transmitted waves allows the reconstruction 
of the scattered wave packets so that the analytical conditions for the stationary phase 
principle applicability can be totally recovered. For such a colliding configuration, the phase 
time gives the exact position of the center of mass of each symmetrically scattered wave 
packet. Consequently, we can have a realistic idea of the magnitude of the dwell time for the 
case in which the energy of the incident particle is smaller than the barrier potential energy 
(tunneling configuration). In spite of the theoretical focus, the results here obtained apply 
only to such configurations which should deserve further attention by experimenters, while 
the existing experiments report the inaccurate results on non-symmetrical configurations. 

In general lines, the SPM allows one to describe the movement of the center of a wave 
packet constructed in terms of a momentum distribution g{k — ko) which has a pronounced 
peak around kg. Assuming that the characteristic phase of the propagation varies smoothly 
around the maximum of g{k — A;,,), the stationary phase condition enables us to calculate 
the position of the peak of the wave packet (highest probability region to find the prop- 
agating particle). With regard to the standard one-way direction wave packet tunneling, 
in which we consider a rectangular potential barrier V{x), V{x) — Vo ii x e [—L/2, L/2] 
and V{x) — ii x [—L/2, -L/2], it is well-known [25] that the transmitted amplitude 



3 



T{k, L) — \T{k, L) \ exp [iQ{k, L)] is written in terms of 

„4 



and 



2k'^ — w"^ 

e{k, L) = arctan \ tanh [2 p{k) L] ) , (2) 



for which we have made exphcit the dependence on the barrier length L, and we have 

1 1 

adopted p{k) = {w'^ — k'^)^ with w = {2mVo)'^ and H = I. The above result is adopted 
for calculating the transit time of a transmitted wave packet when its peak emerges at 
X = L/2, 



tT 



m dQ{k,a) 



2m L ( sinh (a) cosh (a) — [2 — w'^) k^ a 



k=kn 



koa \ 4 fcg {w^ — k^) sinh^ (a) 



(3) 



Here we have defined the parameter a = wL^ (1 — kljw'^) and k^ as the maximum of a 
generic symmetrical momentum distribution g{k — ko) that composes the incident wave 
packet. By following our previous analysis [25], it is well-established that, due to the filter 
effect, the amplitude of the transmitted wave is essentially composed by the plane wave 
components of the front tail of the incoming wave packet which reaches the first barrier 
interface before the peak arrival [27]. We have shown that the cut off of the momentum 
distribution at /c ~ (1 — 6)w increases the amplitude of the tail of the incident wave so that 
it contributes so relevantly as the peak of the incident wave to the final composition of the 
transmitted wave. Independently, due to the novel asymmetric character of the transmitted 
amplitude g{k — ko)\T{k, L)\, an ambiguity in the definition of the arrival/ transmitted time 
is created [25]. In the framework of the multiple peak decomposition [26], we have suggested 
a suitable way for comprehending the conservation of probabilities where the asymmetric 
aspects previously discussed [25] could be totally eliminated. In order to recover the scattered 
momentum distribution symmetry conditions for accurately applying the SPM, we assume 
a symmetrical colliding configuration of two wave packets traveling in opposite directions. 
By considering the same rectangular barrier V{x), we solve the Schroedinger equation for 
a plane wave component of momentum k for two identical wave packets symmetrically 
separated from the origin x — 0. By assuming that (f)^^'^\k,x) are Schroedinger equation 
stationary wave solutions, when the wave packet peaks simultaneously reach the barrier (at 



4 



the mathematically convenient time t — —{mL)/{2ko)) we can write 

4>\^''^{k,x) = exp[±iA;x] + i?^(«)(A;,L)exp[=Fi/cx] x < -L/2{x > L/2), 

(/)^(«)(/c,x) = <J 4>^^''\k,x) = 7^(«)(;fc)exp[^px] +/3^(«)(/e)exp[±px] -L/2<x<L/2, 

4>s^''\k,x) = r^(«)(A;,L)exp[±zA;x] x > L/2{x < -L/2). 

where the upper(lower) sign is related to the index L{R) corresponding to the incidence on 
the left (right )-hand side of the barrier. By assuming the conditions for the continuity of 0^'^ 
and their derivatives at x = —L/2 and x = L/2, after some mathematical manipulations, 
we can easily obtain 

^' ^ ^ \l-exp[2p(A;)L]exp[z2 6'(A;)] J ^' 

and 

r (/c,L)-exp[ ^^^J|i_e3,p[2p(A;)L]exp[i2e(/c)]/' 

where 

^ ' 2k^ -w^ ^ ' 

and R^^'^\k, L) and T'^^^\k, L) are intersecting each other. 

Since the above colliding configuration is spatially symmetric, the symmetry operation 
corresponding to the 1 <-> 2 (or L <-> i?) particle exchange can be parameterized by the 
position coordinate transformation x —x. At the same time, taking advantage from the 
notation that we have adopted, it is easy to observe that 

</.^(«)(A;,x) = ^^,i%{k,x) = ^fi^Uk,-x) = r^'^\k,-x) (7) 

where the L ^ R interchange is explicit. Consequently, in case of analyzing the collision of 
two identical bosons, we have to consider a symmetrized superposition of the L and R wave 
functions, 

(j)+{k,x) = (t>''{k,x) +(j)"-{k,x) = (t)'^{k,-x) + (j)''{k,-x) = (j)+{k,-x). (8) 

Analogously, in case of analyzing the collision of two identical fermions (just taking into 
account the spatial part of the wave function), we have to consider an antisymmetrized 
superposition of the L and R wave functions, 

(t)-{k,x) = 4>^{k,x)-4>''{k,x) = (p''{k,-x)-(p''{k,-x) = -(p_{k,-x). (9) 
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Thus the amphtude of the re-composed transmitted plus reflected waves would be given by 
R^''^{k, L)+T'^'^{k, L) for the symmetrized wave function 0+ and by R^''^{k, L) - T'^'^{k, L) 
for the antisymmetrized wave function Reporting to the previously introduced procedure 
[26] that we call multiple peak decomposition, for such a pictorial symmetrical tunneling 
configuration, we can superpose the amplitudes of the intersecting probability distributions 
before taking their squared modulus in order to obtain 

^ ' ' ^ ' ' ^ \ 1 ± exp [p(A;) L] exp [z 0(A;)] J 

= exp{-z[A;L-(^±(A;,L)]} (10) 

with 

where the plus sign is related to the results obtained for the symmetrized superposition 
and the minus sign is related to the antisymmetrized superposition. Independently of the 
symmetrization characteristic of the wave function, we observe from Eq. (10) that, differently 
from the previous standard tunneling analysis, by adding the intersecting amplitudes at 
each side of the barrier, we keep the original momentum distribution undistorted since 
\R^''^{k,L) ± T'^'^{k,L)\ is equal to one. At this point we recover the most fundamental 
condition for the apphcability of the SPM . It allows us to accurately find the position of the 
peak of the reconstructed wave packet composed by refiected and transmitted superposing 
components. The phase time interpretation can be, in this case, correctly quantified in terms 
of the analysis of the novel phase f^i^k, L). 

Prom this point we have opted for performing the analysis of the symmetrized superpo- 
sition from which, for simplicity, we have suppressed the plus sign from the notation. We 
firstly notice that the novel scattering amplitudes g{k — ko)\R^'^ -\- T^'^| g{k — kg) still 
remain symmetric and, from the SPM, the scattering phase time, results in 

(a) mdip{k,a) 2m L w"^ smh.{a) + ak"^ 

'^''^~~k dk ~ ka 2k^-w'' + w^ cosh (a) ^ ' 

where k /co, with a previously defined [37]. Hereafter, we cannot differentiate the tunnel- 
ing from the reffecting waves for such a scattering configuration. We have simply demon- 
strated that the transmitted and refiected interfering amplitudes results in a unimodular 
function which just modifies the envelop function g{k — k) by an additional phase. The pre- 
viously pointed out incongruities which cause the distortion of the momentum distribution 
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g{k — k) are completely eliminated in this case. The point is that the old phase Q{k,L) 
(Eq. 2) appears when we treat separately the momentum amplitudes T{k,L) and R{k,L), 
which destroys the symmetry of the initial momentum distribution g{k — ko) by the presence 
of the multiplicative term T{k, L) or R{k, L), and the novel phase (p{k, L) appears only when 
we sum the tunnneling/scattering amplitudes so that the symmetrical character of the initial 
momentum distribution is maintained (due to the result of Eq. (10)). It requalifies the SPM 
for accurately computing the time dependence of the position of the peak of a wave packet. 
At the same time, one could argue about the possibility of extending such a result to the 
tunneling process established in a standard way. We should assume that in the region inside 
the potential barrier, the reflected and transmitted amplitudes should be summed before we 
compute the phase changes. Obviously, it would result in the same phase time expression as 
represented by (12). In this case, the assumption of there (not) existing interference between 
the momentum amplitudes of the reflected and transmitted waves at the discontinuity points 
X = —L/2 and x = L/2 is purely arbitrary. Consequently, it is important to reinforce the 
argument that such a possibility of interference leading to different phase times is strictly 
related to the idea of using (or not) the multiple peak (de) composition in the region where 
the potential barrier is localized. To illustrate the difference between the standard tunneling 
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FIG. 1: Normalized times for the standard one-way direction tunneling and the symmetrical scat- 
tering process. These times can be understood as transit times in the units of the classical traversal 
time Tfc = {mL)/k. Both present the same asymptotic behavior. 
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phase time tip^ and the symmetrical scattering phase time ^ we introduce the parameter 

n — k^/w'^ and we define the c/assica/ traversal time = {mL)/k. The normahzed phase 

times can then be written as 

(g) _ 2tu [ cosh (g) sinh {a) - an{2n - l) \ ^(g) _ ^r^ [ na + sinh (a) | 

^ ~ a \ [4n(l-ra) + sinh'(a)] j ^'"^ ~ a \ 2n - 1 + cosh (a) j ^ 

which are plotted in Fig.(l) for some discrete values of n. 

At this point, could one say metaphorically that the identical particles represented by 

both impinging wave packets spend a time equal to tr,^ inside the barrier before retracing its 

steps or tunneling? The answer is in the definition of the dwell time for the same symmetrical 

colliding configuration that we have proposed above. The dwell time is a measure of the time 

spent by a particle in the barrier region regardless of whether it is ultimately transmitted 

or reflected [3] , 

tD = ^ dx\4>,{k,x)\' (14) 

J-L/2 

where jin is the flux of incident particles and 4>2{k,x) is the stationary state wave function 
depending on the colliding conflguration that we are considering (symmetrical or standard) . 
To derive the relation between the dwell time and the phase time, we reproduce the vari- 
ational theorem which yields the sensitivity of the wave function to variations in energy. 
After some elementary manipulations of the Schroedinger equation [28], we can write 

^t^=j_Ad*^_^._^y (15) 

^ ^ 2m dx \dE dx dEdx J ^ ' 



+l/2 

(16) 

-L/2 



Upon integration over the length of the barrier we flnd 

In the barrier limits {x — ±L/2), for the symmetrical conflguration that we have proposed, 
we can use the superposition of the scattered waves to substitute in the right-hand side of 
the above equation, 

4>[^''\k,x) + ct>f''\k,x) 



4>{k,x) 



-L/2{+l/2) ^ 

= exp [zbi fcx] + exp [^i fcx + i ((^(fc, L) — fcL)] (17) 



By evaluating the right-hand side of the Eq. (17), we obtain 

dk dip m /■+^'^^ , I , .,2 Im[ex.-p {iif)] dk 

dEdk=jL,, d.|</>.(fc,x)r h^^- (18) 



The first term of the above equation is the phase time or the aforementioned group delay 
i^^. The second term leads to the exphcit computation of the dwell time. By respecting 
the continuity conditions of the Schroedinger equation solutions, in the barrier region we 
obtain a stationary wave symmetrical in x, 

= \/2(/3 + 7)cosh[p(A;)a;], (19) 
which, from Eq. (14), leads to 

(a) ^ 2 Tfc n sinh (a) + a 

a 2n - 1 + cosh (a) ■ ^ 

The self-interference term which comes from the momentary overlap of incident and reflected 

waves in front of the barrier is given by [38] 

(a) _ /m[exp (zy)] dk _ m sin (y) _ 2Tfc (1 - n) sinh (a) 

k dE~ k^ ~ a 2n-l + cosh (a) ' ^ ' 

In conclusion, the dwell time is obtained from a simple subtraction of the quote self- 
interference delay ti^ip from the phase time that describes the exact position of the peak 
of the scattered wave packets, t^,^ = to,ip + ti^^, as we can notice in Fig. 2. Adopting the 
c/a55«ca/ traversal time r^. = {mL)/k in place of the Earner time r„ = {mL)/w = T^.^Jn 
for normalizing the results illustrated in Fig. 2 may clarify some important aspects. Let us 
assume that the amplitude of each separated transmission coefficient T prevails over the 
amplitude of each refiection coefficient i?, i. e. |T|^ > — > \r\^ > 1/2. For satisfy- 
ing such a requirement the Eq. (1) gives {wL)/{2^Jn) < wLsmh.{a)) /{2^Jn) < 1. At the 
same time, the possibility of accelerated tunneling transitions with respect to the traversal 
c/asszca/ course, i. e. < r^, occurs only when a/2 > (a/2) tanh (a/2) > 1. Since 

a = wLy/1 — n, the intersection of the "weak version" of both of the above constraints, 
{wL)/{2^Jn) < 1 and a/2 > 1, leads to n > 2, which definitely does not correspond to a 
tunneling configuration. This analysis leads to a very important conclusion: in the region 
where the one-way direction transmission coefficient prevails over the refiection coefficient, 
the wave packet transmitted components should tunnel with a retarded velocity with respect 
to the classical velocity since, in this case, > r^.. The supposed accelerated transit of 
the tunneling wave packet, and therefore, superluminal velocities and the Hartman effect, 
would never occur for I TP > 
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FIG. 2: Exact phase time (solid line), self-interference delay (dotted line), and the dwell time 
(dashed line) as a function of the normalized energy n = k'^/w^ oc Eq/Vq for the symmetrical (black 
line) and the standard one-way direction (red line) configuration for wave packets colliding with 
a rectangular potential barrier. These times are normalized by the barrier time = {mL)/w in 
the first plot and by the classical traversal time = {mL)/k in the second plot. Here we have 
adopted wL = Att for a = wL\/\ — n. 

To summarize, we have reported about a way of comprehending the conservation of prob- 
abilities [25, 26] for a very particular tunneling configuration where the asymmetry and the 
distortion aspects presented in the standard case were all eliminated. As a result the phase 
time could be accurately calculated in order to give the exact position of the scattered wave 
packets (or particles), for which, in the standard case, the position of the peak is shifted. 
Essentially, we have claimed for relevance of the use of the multiple peak decomposition 
technique in obtaining the exact relation between two distinct scattering/tunneling time 
definitions: the phase time and the dwell time. In spite of quoting the superluminal inter- 
pretation, our discussion concerned with the definition of the strict mathematical conditions 
stringing the applicability of the stationary phase principle in deriving transit times. Even 
with the introduced modifications, our results corroborate with the analysis [11, 21] that 
gives an answer to the paradox of the Hartman interpretation [22]. In certain sense, we 
accept the fact that one should keep in mind that the Hartman Effect, even in its more so- 
phisticate consequences) appears to have been experimentally verified [29], in particular, for 
opaque barriers and nonresonant tunneling [30] and also reproduced by numerical simula- 
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tions and constrained theoretical analysis [31, 32]. In general lines, there have also been some 
trying of yielding complex time delays for tunneling analysis, ultimately due to a complex 
propagation constant where the supposition of superluminal features is considered artificial 
since the transmitted peak is not causally related to the corresponding incident peak. This 
semi-classical method makes use of complex trajectories which, in its turn, enables the defi- 
nition of real traversal times in the complexified phase space [7, 33, 34] . Yet concerning the 
subsequent theoretical perspectives in the above frameworks, the symmetrical colliding con- 
figuration also offers the possibility of exploring some problems involving soliton structures. 
In particular, all the above arguments suggest that the idea of complexifying time should 
be investigated for some other scattering configurations, which reinforces the more general 
assertion that the investigation of wave propagation across a tunnel barrier has always been 
an intriguing subject which is wide open both from a theoretical and an experimental point 
of view. 
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